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A THEORETICAL STUDY OF THE MOMENT ON A BODY IN A COMPRESSD3LE FLUID 

By Carl Kaplan 



SUMMARY 

The extension to a compressible fluid of LagaUy's 
theorem on the moment on a body in an incompressible 
fluid and Poggi's method of treating the flow of compres- 
sible fluids are employed for the determination of the 
effect of compressibility on the moment on an arbitrary 
body. Only the case of the two-dimensional subsonic flow 
of an ideal compressible fluid is considered. 

As examples of the application of the general theory, 
two well-known systems of profiles are treated; namely, 
the elliptic profile and the symmetrical Joukowski profiles 
with sharp trailing edges. 

The effect of compressibility on the position of the 
center of pressure is also discussed. In order to deter- 
mine this effect, if is necessary to calculate the additional 
circulation induced by the compressibility of the fluid for 
both the elliptic and the Joukowski profiles. For these 
two types of profile, the centers of pressure in the com- 
pressible and the incompressible fluid* are found to coin- 
cide for a definite and fairly small angle of attack, which 
is essentially dependent on the thickness coefficients. For 
angles of attack less than this neutral angle, the center of 
pressure in the compressible fluid is farther from the nose 
and, for angles of attack greater than the neutral angle, 
nearer to the nose than the center of pressure in the in- 
compressible fluid. 

Several numerical examples of both the elliptic and the 
Joukowski profiles are given. The results show that, 
although the effect of compressibility on the moment and 
on the lift may be large, the effect on the center of pressure 
for conventional profiles is negligible. Thus, for a 
Joukowski profile, the maximum thickness of which is 
equal to IS percent of the chord, the center of pressure 
moves toward the nose a distance equal to only 0.19 per- 
cent of the chord, where the angle of attack is 6° and 
ojc 0 =0.70. 

NOTATION 

GENERAL SYMBOLS 

§, tj, rectangular Cartesian coordinates in the 

plane of the obstacle. 
x, y, rectangular Cartesian coordinates in the 
plane of the circle. 
r=£+£ij t z=x+iy 



r, 6, polar coordinates in the plane of the circle. 

B, radius vector of a point far removed from 
the obstacle and also of a point far re- 
moved from the corresponding circle. 

r 0j radius of circle into which the profile is 
mapped. 

r 

v t , Vs, components of the velocity in the radial 
and the circumferential directions in the 
z plane. 

v, magnitude of the velocity in the plane of 

the obstacle, 
c, magnitude of the local velocity of sound. 
p, density of the fluid. 
p, static pressure of the fluid. 
% 4>, poi Pa, corresponding magnitudes in the undis- 
turbed part of the fluid. 
y, ratio of the specific heats (c p fc,). 

M=Q^ » square of the Mach number. 

jS, angle of attack. 

r, circulation about the obstacle. 



K= 



L=p {i v 0 T, 
M, 



lifting force on the obstacle, 
moment on the obstacle. 
c e , c u centers of pressure, respectively, in the 
compressible and the incompressible 
fluids. 

Subscripts c and i refer, respectively, to 
the compressible and the incompressible 
fluids. 

Av a , Avs, additional velocity components due to 
compressibility at a point P(R, 5) far 
removed from the circle in the z plane, 
complex velocity potential of the fluid. 



w, 
dw 

ap 

A a ,A a , 



complex velocity in the plane of the ob- 
stacle. 

complex and conjugate complex coeffi- 
cients, respectively, in the power series 

development of (see equation (12)). 
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SYMBOLS PERTAINING TO ELLIPTIC PROFILE 

c, semifocal distance. 

t, thickness ratio of ellipse (ratio of semi- 
minor and semimajor axes). 

tank a= t 

M, N, functions of the thickness coefficient t only 
(see equation (24)). 

SYMBOLS PERTAINING TO JOUKOWSKI PROFILES 

e, thickness coefficient (see fig. 4). 

h ~ l+€ 

l+« 

/, J, (see equation (26)). 
M, N, functions of the thickness coefficient e only 
(see equation (29)). 

DERIVATION OF THE FORMULA FOR THE MOMENT 

Theodorsen's method (reference 1) of extending 
Lagally's formula for the force on a body in an incom- 
pressible fluid to a compressible fluid may be used to 
obtain the corresponding formula for the moment. 
The body is fixed in an infinite two-dimensional stream 
of a Motionless compressible fluid flowing uniformly in 



V 




i 



Figure 1.— Directions of tangent and normal on a contour. 



the far field. Then the moment exerted on the body 
with respect to the origin of the coordinate system is 
given by (fig. 1) 

M=—f t [r n]p ds (1) 

where the positive direction of the unit normal vector n 
is taken from the boundary into the fluid, p is the 
pressure of the fluid, [r n] is the vector product of the 
radius vector r and the unit normal vector n, and ds is 



the element of length along the profile taken positively 
in the direction of the unit tangent vector t. With the 
directions of the unit tangent vector t and the unit 
normal vector n so chosen that they form a right-hand 
system, a positive value for the moment corresponds to 
a counterclockwise rotation. 

According to one of the generalizations from Gauss' 
theorem, it may be shown that 

f, curl ptdS=f i [v n]p ds-\-f„[r n]p ds 

where dS is the element of surface in a region S in- 
cluded between the obstacle i and an arbitrary curve 
o enclosing it. 
Now 

curl 2>r=[grad p i] 

and, from the Euler equations of motion for steady 
flow (reference 1), 

grad p= — 1/2 p grad »* 

Therefore 

curl pr= — 1/2 p [grad v* r] 

Equation (1) then becomes 

M=f e [r n] p tfe+1/2/, p [grad tfi r] dS 

Since the outer boundary is arbitrary, it may bo chosen 
to be a large circle so that the vector product [r n]=0. 
Hence 

il/=l/2/ aP [grad^r](iS 

But 

grad v*=2 (v grad) v+2 [v curl v] 

and 

[(v grad) v r]= (v grad) [v r] 

Therefore 

M-f. p (t grad)[v r]dS+f. P {[v curl v]r}dS (2) 
According to Gauss' theorem, if F denotes a scalar 
function, then 

f„ div P FydS=-fi P F(v n)ds-f 0 P F(v n)ds 

But 

div P F\=F div pv+p(v gr&d)F 

Hence 

f, p(v grad)i%S= — f, F div P vdS 
-f, P F(v n)ds-f B P F(v n)ds 
From the manner in which F occurs in this formula, 
it is clear that the formula remains valid if F is re- 
placed by a vector point function, say [v r]. That is, 

f* p(t grad) [v r]dS= f, [r v] div pvdS 

+ft p[r t](v n)ds+f 0 P [r v](v n)ds 
Substituting from this equation into equation (2), it- 
follows that 

M= fo p [t v] (v n)ds+ f, [r v] div pvdS (3) 

- Jl P {r[vcmi \]}dS 
where y* f p[r v](v n)as=0 since (v n)=0 everywhere on 
the obstacle. 1 

In the problem considered herein, the compressible 
fluid is of uniform velocity c 0 in the undisturbed stream 

i This formula for the moment may be extended to three dimensions by letting. 
t denote the region of Sow between the body and a large sphere enclosing It; it Is 
then replaced by do, the element of surface on the outer sphere, and <f S Is replaced 
by &t, the element of volume in the region of flow r. Also, the moment now has 
three components and must be written as a vector, I. e., M. 
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and, furthermore, is assumed to be irrotational and 
free of sources in the region of flow S. The condition 
for irrotational motion is simply that curl v=0, and 
the absence of sources means that div pv=0. There- 
fore, equation (3) becomes 

M=f eP [tY] (vn)<fc (4) 

where it is recalled that the outer boundary has been 
assumed to be a large circle of radius R. The vector 
product [r v] is then equal to Rv v and the scalar prod- 
uct (v n) is equal to — v B since n= — E/H. Therefore 

M=-R 2 f, 2T pv B v r d<p (5) 

This simple expression for the moment about the origin, 
of the resultant force acting on the body, can be directly 
obtained by considering the rate at which angular mo- 
mentum passes out of the region included between the 
outer circle and the body. Thus, the momentum per 
unit time passing normally through the element ds is 
p Rv e B d<p and the arm is R sin vR. The angular mo- 
mentum is therefore pR 1 v sin vR v B d<p or pR 2 v B v r dp 
and equation (5) follows. 

GENERAL DEVELOPMENTS 

It may be well to emphasize at this point that the 
main problem of this paper lies in obtaining useful 
expansions for the velocity components r B and v f of a 
compressible fluid. 

These expansions are obtained by making use of 
Poggi's conception of compressible flow. Thus the 
basic differential equation for the steady flow of a 
compressible fluid may be written as (reference 2) 



The expression on the left-hand side is div v, so that the 
expression on the right-hand side may be considered to 
represent a source distribution of a strength given by 



4tc i ( 



In the plane of the circle into which the profile of the 
obstacle is mapped by a suitable conformal transforma- 
tion, the strength of the source distribution may be 
written as 



where 



r, 6 are the polar coordinates of a point in the 
plane z(=x+iy) of the circle. 

X=— , where r 0 is the radius of the circle into 
r 

which the profile is mapped. 
5r =_&* P(== _I|* where * is the velocity 
potential of the flow. 



v, the magnitude of the velocity of the fluid in 

the plane of the profile. 
c, the magnitude of the local velocity of sound. 

Poggi's method of approximating the flow of a per- 
fect compressible fluid past an obstacle is based on the 
fact that the incompressible fluid is a good first ap- 
proximation. A first approximation for the source 
distribution given by equation (6) is therefore obtained 
by substituting for v r> u, and & the values belonging 
to the incompressible fluid. It is further assumed that 
the Mach number r 0 /c 0 is small so that only terms 
involving the lowest order of r 0 /c 0 are to be considered. 
This assumption then limits the application of the 
analysis to stream velocities small in comparison with 
the velocity of sound. The disturbances to the main 
flow due to the presence of a body in the fluid may, 
however, be large. 1 

The density p and the pressure p of the fluid are 
then determined by the following equations (reference 
2): 



=A [l + l/2p(l-g) + . . .] 
p=constant-|p 0 p (! 2 -|jt:po»f ii +|j^r--|- 



(7) 



where the adiabatic equation of state has been adopted. 
Also, p=(Pb/c 0 )* and v f and v, are, respectively, the 
velocities in the incompressible and the compressible 
fluids. 

It follows from equation (5), with the stipulation that 
only terms mvolving the square of the Mach number are 
to be retained, that: 



(8) 



It is clear that the evaluation of the second integral 
presents no difficulties, for it involves only expressions 
of the well-known incompressible fluid-velocity com- 
ponents. The difficulty of the problem lies mainly in 
determining the velocity components v B and v v for a 
compressible fluid. Since the first integral is taken 
around a circle whose radius R may be infinite, it is 
obvious that the developments for v B fv 0 and v v fv 0 in the 
neighborhood of infini ty need not go beyond the l/i2 2 
terms. 

In order to obtain the series for Vg/vo and v v fvo, it is 
expedient to consider first the effect of a single source 
of unit strength situated at a point Q(r t 6) in the plane 
of a circle. In the presence of a circular boundary of 
radius r 0 , the velocity induced by a unit source at any 

' In the approximation made by Glauert (reference 3), n Is not small compared 
with ri hut the disturbances to the main flow tt are assumed to be small. 
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point P(R, 5) external to the boundary is given by 
(fig. 2): 




FiouM 2— Image of a simple source with regard to a olrcle. 

dw /_1 . 1 J_\ _ w + . v 

d~2 \2 P —2<}2p—2a 2 P ) " 

where £ is the point inverse to Q with respect to the 
circle; 



z P =Re iS ; ZQ=re {) ; and z 8 =-je ti . 



Then 



v x —w v = 



1 



1 



or suace 



for r 0 <r<B 



* Re«-re^ Re(i __rl eit Re« 

T 



(— »fi+i«j)e -<s = — v z +iv v 



1 

Vr =2R 



I — 



X 2 



1-2-^ cos #-5)+^ 



l-Xp 2 X* 



+ 2~2? 1-2X P X cos (0-5) + Xp J X J 



sin (0-5) 



1-2%- cos (0-5) 



1 X P X sin (0— 5) 

~R l-2XpX cos (0-5) + Xp*X s 



and for i?<r< 0= 



v a =- 



2i? l-2^- cos (0-5) 

Ap Ap 



1 l-Xp 2 X 2 
+ 25l-2XpX cos (ff-SJ + VX*" 



Pj= — 



sin (0-5) 

Ap 



R l-2^- cos (0-5)+S 

Ap Ap 

__1 XpX sin (0—5) 

R l-2XpX cos (O-Sj-l-Xp'X 1 

where X=r 0 /r and \ P —r 0 /R. 
Then, making use of the expansion 

T— s — ia , »=l+2Sg" cos n(0-5) 
1— 2g cos (0— 5)+^ iTTi 

it follows that the components of the velocity induced 
at any point P(J?, 5) by the source distribution given 
by equation (6) are: 



+ 



and 



Jf* sin (0-5) 



x8 

*o 4xJx, Jo ! _ 2 X^p cog ( ^_ 5) + X^ 
^"4ttJo Jo j 



X v " ' X> 

sin (0—5) 
.2 A cos 0-J)+£ 

Ap Ap 



cfX c?0 



+ 



47rJ 0> /o 1— 2 XpX cos (0—5) + Xp 4 ^ v ' 



where 



\Vo d\ «o X £»0 / f 



The next step is to obtain the Fourier series expansion 
of the function F for an arbitrary profile. The profile 
is derived from a circle by means of the general con- 
formal transformation 

where. the coefficients a t are, in general, complex. In 
order to obtain the velocity of the fluid in the plane of 
the profile, it is first necessary to determine the velocity 
potential of the flow past a circle (fig. 3). Thus, sup- 
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pose that the undisturbed flow is of velocity ro inclined 
at an angle /S to the negative direction of the real axis 
and that the circulation is denoted by T. In terms of 



Q(n,e) 




> x- 



Figure 3.— Eelatlon between r planes employed In general confonnal transformntton. 

the complex coordinate z, with origin at the center of the 
■circle, the potential function of the flow past the circle is: 

Also, the variables 2 and z' are connected by the equa- 
tion 

s'=z+ me 1 ' 

Therefore 

dw = dw dz dz' W 2* 4 
2T dz dz'~dj <ii _2<Z2 

1 i 74 i 7 ^ 



2x2 



Substituting for z' in terms of 2 and expanding in 
descending powers of 2, this expression, becomes: 

if= e ' fl +2S 0 i 2+^-^4+--- <«> 

It follows then, that: 

«oC?f £o2 n 

where the A*'s are, in general, complex. 

In order to introduce the angle 6— S, the general 
term of this series is multiplied by e~ ini je~ tne . 
Then 

and 



1 dw_^ AfiW 



The real and the imaginary parts of A t Aje l ^~ m are, 
respectively: 



and 



Let 



A t A J e t ti-*> l +A t A J r tt t- k,t 
2 



AtA/v-^-AiA/i-w-™ 



j—k=n 

where n changes sign when j and 1c are interchanged 

and therefore takes values from — » to ». 

Then 

and the expression for c 2 /^ 8 , a real quantity, becomes: 
^=2 S r^»cosn(5-«)+5,. t sinn(0-5)]X» +2i (13) 
where 

2r 0 n+i * 



and 



■n .AtA^'-AiA^-™ 



From the definition of 71, it follows that: 

A—n,t = An,k 

and 

The terms of equation (13) can therefore be grouped in 
pairs, thus: 

~=f:\ u A 0 , t +2^ cos n(ff-5) fj A.*X* +2 * 

V 0 k^o n-1 i^O 



+ 2Z) sin n(fl-5)2-B»^X» +M 

n=l t^O 

Also, from equation (11), 



(14) 



and 



where 



- r =-(l— X s ) cos (0-8) cos («+jS) 
+ (1-X i! ) sin (0-5) sin (5+j3) 



-=(1 + X ! ) cos (0-5) sin (5+/3) 
+ (1 + X 2 ) sin (0-5) cos (5-f-jS) + XK: 

r 



(15) 
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Then, by means of equations (14) and (15), it becomes 
a simple matter to obtain the Fourier series for the 
function F. When the Fourier series for F is substi- 
tuted into equations (9) and (10) and all terms con- 
taining powers of Xj» higher than the second are neg- 
lected, it follows, after some rather tedious but element- 
ary integrations, that: 

^*=-£\ P [Aw cos (5+j3)-F w sin (S+/3)] 
-ftA,! cos (5+ p) +KB lfi W log Xp 
-MXpjj^ cos (5+/3)+^t 2 ,o cos (5+0) 

-^ S in( 5+ «] + fv|;|^ B , a 



+|x^cos(8+«Sj±^A» +1 



1 n-l 



A. cos (8+0)-B 1§tt sin (5+g) 
n+1 



(16) 



and 

At), 



X^.0 sin (S+P)-B lfi cos (5+ (3)1 



-gWw sin (J+ft+KAJV log Xp 



+ M X i , 2 [^, 1 sin («+/3)+§Aosin (8+|3)+^x.o] 
+|V-(^)gJ^, +1 



+fx^S 



A.» sin (8+|3)+-g».»cos (5+j3) 

71 + 1 



(17) 



Use has been made of the following definite integrals 
in obtaining Av s (v 0 : 

f 2 * sin (3-5) sin n(d-5)de TO if n=0 
Jo 1-2A cos (0-5)+ A 2 [,rhr- l Un^l 

p* sin (g— 5) coa.7i(g— 5)dfl _ n 
Jo l-2Acos (0-5) +A 2 

The velocity of the compressible fluid in the plane of 
the circle then has the following components: 



and 



(18) 



The calculation of the moment is facilitated by 
expressing the velocity of the compressible fluid in 



complex form. Thus, as in the case of an incompres- 
sible fluid, a complex velocity is defined in the following 
way: 

v<\dz/ e \ V 0 Do A 

Then equation (8) for the moment assumes the Blasius 
form: 

M.~i*,.R.P.$(i£)>, 

>«' P 'f('-S(«)> ™ 

where the integrals are taken around a circle whose 
radius R approaches infinity. 

Now,' from equations (13), (16), (17), and (IS), the 
expression for l/» 0 (dw/dz) c may be written as: 



Hp) =e 



iKr 0 r 0 *e-*<> 

' Z 3 s " 



^A^A^+A\A^\ 



+2(A 1 X l e-*-tKr 0 Z ft A 1 )p log Jf+^e* 

. .„! » 2n+l M,ti 
2e-<^ n+1 A^An+r , 2e^ 1 A^j ] 



z* £in{n+2) r 0 *» 
Also, from equation (12) for ^ 

_ A=e^, Ai=iKr 0 , A^bV^-rfe-* 
where ai=b 2 ^ 
Therefore 



i 5 log B + ^ 



IF 



" 2» 22 s ^ -* 



, 6V^>-r 0 2 e-^ , iE^ 2n+l 
+ <J 3 s ^z'&inin+l) ro 2 - 1 

z* t=in{n+2) r 0 tn + s 2 f=fc+l 



A THE ORETICAL STUDY OF THE MOMENT ON A BODY IN A COMPRESSIBLE FLUID 



587 



It may be noted that this expression is not an ana- 
lytic function of 2, since powers of B occur in some of 
the terms. 

The first integral on the right-hand side of equation 
(19) for M e then becomes: 



~fW R. P 



1 . dw\ 2 dz 



= - lA3 ^R.P.^-^(l+— + 



z 3 



=2trb 1 p 0 v a 1 sin 20S+y)+poi' 0 rm cos 03+?) 
+|[3 1 r6 ! p 0 i;o ! sin2()3+r) 

+^r i E.P. ^£+} ) ^ 



-2ta,Co 2 R. P. iV*2> 



a+2 



0 



if=l(a+l)r a 2 » 
and the second integral: 

=|[2ir6*poi?<,* sin 2(/3+Y)+PofoTm cos QS+jO] 
Therefore 

M e =M i (l+f)+PoCbAr7n cos (0+*) 

+^r,B.P. ^+ 1 1 ) ^ 



-2tpo« 0 * R. P. is 2 *" 



^(7i+l)r 0 2 -»J 



(20) 



where AI t =2icb 2 p a v Q i sin 2(^+7)+Po«or { m cos OS+v) is 
the moment due to an incompressible fluid. The cir- 
culation T consists of the incompressible part r ( and 
the additional circulation Ar related to the additional 
flow. 

Equation (20) for the moment M e is applicable to 
arbitrary profiles, but it contains two infinite series 
that cannot, in general, be replaced by closed forms. 
In the case of a profile for which the conformal trans- 
formation to a circle contains a finite number of terms, 
the two infinite series may, however, sometimes be 
replaced by elementary functions. In the following, 
two such systems of profiles will be discussed, the 
elliptic and the JoukoWski profiles. 

THE ELLIPTIC PROFILE 
It is well known that the Joukowski transformation 



maps the circle of radius c/2 with its center at the 
origin of the z plane into a line segment (— c, 0; c, 0) in 
the f plane. Also, circles concentric with the base 
chtole are transformed into a family of confocal ellipses 
with common focii at (— c, 0) and (e, 0). If r 0 (>c/2) 
denotes the radius of one of these circles, then the 
semimajor and the semiminor axes of the ellipse into 
which it is transformed are, respectively, r 0 -t-c*/4ro and 
r„— c*/4r 0 . The thickness ratio t is then defined as 



4r 0 l-ir 2 



r °+47 0 



or 



where 



1-t 

1 + t 



2r a 



Now, for the case of the elliptic cylinder, equation (12) 
becomes: 

l rfw = ^2t% z 2 s 
v a dj 



\ + 2*t> 0 z z> J 2 -'^) ~ 2 ~ t z« 



or 



where 



Then 



^ 2rc^ 0 



riB ^ 2n+l „A 4n+l A^A^ +l 

t=in{n+l) ro 2 " _e ^2n(2«+l) r 0 te 



4n+3 



2»+l-£»-2il+g 



Z=l&n+1) (2*1+2) r 0 ta+J 
" / 1 



1 V- s 

lw\2n ' 2n+lJ 
v.2n+1^2n+2/ 



and 



"R P f»T' 27l4"l -^-n^-fl 



=i£> 0 sm 2^S— log j^p-pj 



Similarly 



-^■ia-l-^aic+l 



^fa+lK 2 " ~ ^=i(27i+l)r 0 ta ^ t=i 2nr 0 te - J 

,to-2 "o^"- 



=ie i " , r 0 2 (l-2 ff s cos2/3+ ff «)S2 W+ j 



-i 2n 
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and 

R. P. ie^S- 



n+i 



= -r 0 2 sin 2/5(1-20* cos 2/3+ o*)(Jp, log fizfl-Ji) 



-if ( % 2 sin2i3S- 



Therefore, since m=0 and p=0. for an elliptic profile, 
equation (20) becomes: .. . 

M e =M t \l +|+H(i_ 2ff * cos 2/J+ o(^3 log x^-Ji) 

+K^- log S-2^-^^ (i -4 (21) 

where Mi—^p^c 2 sin 2/3. 

For thick profiles, or for small values of a, this formula 
may be expressed as a power series in <r. Thus, by 
making use of the expansions 

log f±£-2(^+|^+ • • •) 

and 

log (l-<^) = -(^+^+^+ • • •) 
it follows, from equation (21), that 

M e =M[l +g+|(l-2<r a cos 2/3+«*)(|+io*+|^ 

+ • • ■)+KI+i^+f 5 ^+ 2 l ffls + ■ • •)] 

It is seen from this equation that the value of ^'jj^* 

for the limiting case of a circle, for which < = 1 or c= 0, is: 

M.-M t _n 4+5g t a 
3 



M t 2 

On the other hand, equation (21) shows that, for the 
limiting case of a straight-line segment for which f=0 
or <r=l and the angle of attack /3 is finite, 

Mc-Mt 
M t =a 

It is to be noted that, although M ( , and therefore M c , is 
zero when the angle of attack /S is zero, the ratio 
M-M t 



possesses the limiting value 



where K t =2 sin /3, according to the condition that the 
rear stagnation point occurs at the end of the major 
axis of the ellipse. For the two limiting cases of the 
circle and the straight-line segment, this ratio becomes, 
respectively: 

( M-M t \ _4 n 

and 



M c —M t \ ]x 
j/s-o=2 



(23) 



In appendix A, an independent derivation of equa- 
tion (21), with Ki=0, is obtained by a direct integration 
of the pressures over the surface of an elliptic cylinder. 
It shows clearly the superiority of the present general 
method, in which all integrations are performed along a 
circle at infinity. 

It is of interest to investigate the effect of compressi- 
bility on the center of pressure. If the rear stagnation 
point is supposed to be at the end of the major axis, it 
is immediately seen from the second of equations (18) 
that, for 5=ir and \ P =1, 

Since K { =2 sin p, it follows that the additional circula- 
tion AK is given by 



AK- 



The problem of determining the additional circulation 
thus amounts to finding an expression for Avi/vo at the 
stagnation point 8=ir. This calculation is given in 
appendix B and it is shown there that 



AT 
ft 



AK 



(24) 



where M and N are functions of the thickness coeffi- 
cient t only. Table I presents values of M and JV for 
various values of t. 

TABLE I 



t 


M 




0 


1.00000 




.1 


1. 09808 


3.44779 


.2 


1. 19226 


2.48068 


.3 


1.28287 


1.98162 


.4 


L 87929 


1.63730 


.8 


L 45392 


1.39S13 


.8 


1.83803 


1.20008 


.7 


1. 61342 


1.03764 


.8 


L 68916 


.89791 


.9 


1.78188 


.77605 


1.0 


1.88333 


.68667 



For the straight-line segment, Ar/I\ equals infinity if 
the value of /S is other than zero but, for very thin 
profiles and vanishingly small angles of attack, 

r f +Ar n 

Tt i "2 

This result agrees with that given by Glauert (reference 
3), namely: 

Tt±AT_ 1 



r, VT 1 



i+S+ 



In reference 3, Glauert has also shown that the lifting 
force on a body in a compressible fluid is given by 



or 



L e =poVoT 



(25) 
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If c„ and Cf denote, respectively, the centers of pressure 
in the compressible and the incompressible fluids, then 
according to equations (23) and (25) 

That is, for a very thin profile and for vanishingly small 
angles of attack, the center of pressure is unaffected 
by the compressibility of the fluid. For the general 
elliptic profiles, it follows from equations (21) and (24) 
that 



or the motion of the center of pressure away from the 
origin is 0.15 percent of the chord. Again, if the angle 
of attack is vanishingly small, then cos jSs^l and sin /3sij9, 
and 



^=1.14486, -^=1.18174 



Then 



c e —c t _ 



2a 



-0.0039 



c _ e= i+fji+ a-^ «» gg+^Gg lQ g ct4)+ 4 Hffi Ip g T^-2T*-h *°g a-^]} 



l+f(M+i\Tsin s /3) 



It may be shown from this expression that, for any given 
ellipse, there exists an angle of attack /3, independent 
of the stream velocity r 0 , to the first order of for 
which the ratio c e /c t equals unity. Furthermore, if the 

angle of attaok is|^ ater Jthan this neutral value of jS, 
then c/ct is|^^" ter Jthan unity and the center of pres- 
sure c, in the compressible fluid is{^„^ er /f om lthe 
c c [nearer to J 

origin than the center of pressure c t in the incompres- 
sible fluid. Table II presents the neutral values of jS for 
the entire range of ellipses including the straight-line 
profile and the circle. 

TABLE II 



t 


cos 20 


Pneutral 


Degrees 


Minutes 


0 


L00O0Q 


0 


0 


.1 


.9R302 


5 


17 


.2 


.08199 


7 


M 


.3 


.91116 


9 


£3 


.4 


.92173 


11 


26 


.5 


.9034S 


12 


40 


.1 


.S8SS7 


13 


46 


.7 


.87148 


14 


42 


.8 


.869*15 


IS 


22 


.9 


.84471 


16 


n 


1.0 


.£3333 


16 


47 



As a numerical example, consider an elliptic cylinder 
with a thickness coefficient t=H and with 5 0 /c 0 =K- 
In this case, the neutral angle of attack /3 is given by 
cos 2/3=0.90348, or 0=12°4O'. If the angle of attack 
is increased to 15°, say, theD 



and 
Now 

Therefore 



4^=1.20776, ^=1.19345 
Alt L t 



- c =1.012 

Ct 



ci=4(i-f) 



2a 



=0.0015 



or the motion of the center of pressure toward the 
origin is 0.39 percent of the chord. In general, when 
the angle of attack 0 is vanishingly small, the center 
of pressure c c in the compressible fluid is nearer to the 
origin than the center of pressure c { in the incompres- 



sible fluid. Table III gives values for (%~) 
various values of t and for the critical values of p.* 
TABLE in 



for 



t 






( *2« On 


0 


1.00000 


1.00000 


0 


.1 


.73445 


.97960 


-.00459 


.2 


.67808 


.9098* 


-.00602 


.3 


.46077 


.96459 


-.00620 


.4 


.39614 


.96155 


—.00577 


.6 


.83293 


.96038 


-.00496 


.6 


.28944 


.95971 


-.00403 


.7 


.25120 


.96003 


-.00300 


.8 


.22212 


.90087 


-.00198 


.9 


.19714 


.961 IS 


-.00097 


1.0 


.17010 


.06204 


0 



THE SYMMETRICAL JOUKOWSKI PROFILE 

The Joukowski profiles are derived by means of the 
conformal transformation 

which maps the circle of radius a with its center at the 
origin of the z' plane into a line segment (—2a, 0; 2a, 0) 
in the f plane. Any other circle of radius r a (^>a) with 
its center lying on the real axis at a distance ae from 0 
(see fig. 4) and touching the base circle at (—a, 0) is 
mapped into a symmetrical Joukowski profile with a 
sharp trailing edge. The variables z and z r are con- 
nected by the equation 

z'=ea+z 

The circulation T t is chosen in accordance with the 
Kutta-Joukowski hypothesis that the rear stagnation 
point of the flow occurs at the point A of the circle. In 



'The critical value of p Is defined as that value of p at which the maximum velocity 
of the fluid at the surface of the obstacle just attains the local velocity of sound. A 
Est of such values of gfor a set of elliptic cylinders with the angle of attack equal to 
isro is given In table TV of inference 47. 
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where 



terms of the complex coordinate z, the potential func- 
tion of the flow past the circle is 

The complex velocity at the point A is then given by 



2 urn 



Therefore 



r ( =4?rtW) sin 0 qtJ£j=2 sin /9 
y'y 



/==K ' l 2sin2|S^7i(ji+l) 7?^ 



J=-R. PS 



and 



"sin 2j3^(n+ 1)^ 1 
M t =2Tp*>W [1+6(1+6)1 sin 20 



— A 




Figure 4.— Transformation of a symmetrical Jookowakl profile Into a circle. 

The complex velocity in the J- plane is given by 



dw dwdz dz' (g+ro)(ge < g-roe-^)^ ,^ 
dj- Tzdl' 'dj ==v ° z*(.z'+a){z'-a) 
Then, since 

3 '-j-a=z+a(l+«)=z+r 0 and z'—a=z—a(l—e) 
it follows'that 



,* (i + l)'(^-v») 

— axn ; ■ ■** ,- r 



v 0 d}~* ' a(l — t) 

Putting YZf e ~^ an( * Xip^ = ^' 

V 0 dt ^=0Z" £T\ZJ n-0 \Z/ 

Therefore 

i4o=«*, Ai=2tr 6 sin 0, 4j=a l «*— r«V*, 

and 

A=sV" 4 ^(l;e i|, -e" < ') for "53 
Equation (20) giving the moment M„ then may be 
written as 

M-M t _ 6(1+6) AT 

~~ ^ 1+6(1+6) r, 
M r 3 6(i+e) fJ+J /l (26) 



In appendix C, closed expressions are derived for I 
and J. They may be written as follows: 



(l+e)»(I+J)=^^(l-2fc cos 2J3+P) 
[^lo gT ^p-(l+2i+|p+^+^+^)] 
+i(3+2fc)(l-2fc cos 2(3+i 2 )+^(9+5/:)(l-cos 2/S) 



+^(1-^+^(15+54+4^(1 



(27) 



For small values of k, or thick profiles, 

(l+6) l (/+J r )=^J^- 2 (l-2A: cos 20+F)(^+ Jfc 

+ i(3+2*)(l-2fr cos 2j3+P)+i(9+5A-)(l-cos 20) 
Ar)+-^(15+5A:+4A: a )(l-A:) a 



(28) 



It remains to evaluate the ratio of the additional cir- 
culation Ar and the circulation T< of the incompres- 
sible fluid. This ratio may be obtained from equation 
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(22) of reference 2 by evaluating Avifvz for 5=ir and is 
written as follows (see appendix D): 



r, 



2 sin 



hn&),.rl( M+N ^ (29) 



where M and N depend only on the thickness coefficient 
e of the profile. Table IV lists values of M and N for 
several values of e. 

TABLE IV 



e 


M 


N 


0 


1.00000 




.03 


1.03242 


3.76210 


.OS 


1. 06496 


3.78831 


.10 


L 11153 


2.SS5S5 


.20 


1. 21659 


2.17368 


.30 


1.32817 


1.80592 


.40 


1.40925 


1.58348 


.GO 


L 49044 


L 38092 



It is interesting to note that M= 1 for a straight-line 
profile but that N is infinite on account of a term con- 
taining log e. For very thin profiles and vanishingly 
small angles of attack, however, 



r t +Ar 



(30) 



This result agrees with Glauert's well-known formula 
(see section on the elliptic cylinder) 

r t +Ar _ i 
r t VT=^ 5 ' ' ' 

It follows from equations (26) and (27) that, for an 
infinitely thin profile with an angle of attack so small 
that sin jS^S and cos /Ssl, 

With ie=Poffo(r<+Ar) and L^/WaTt, it therefore 
follows from equations (30) and (31) that 

M e _Mt 

or 
c 

■7=1 (as in the case of the ellipse). In general, ac- 
cording to equations (26) and (29), 



In order to gain some idea as to the order of magnitude 
of the movement of the center of pressure due to com- 
pressibility, consider the case of a thin profile, say 
«=0.05. Then equation (32) becomes 

1 + £(7.00142— 5.95991 cos 2/3) 

Ct l+|(2.92162-1.86666 cos2/3) 

Letting c e and c t coincide, i. e., Cc/c t = 1, yields 
cos 20=0.99671 

or a neutral angle 

0=2°2O' 

If the angle of attack is taken to be 0°, then 

c e _ 1+0.52075/1 
c t 1+0.52748 fi 

which shows that c e is nearer the origin than c t . It is 
seen that the movement of the center of pressure is 
very small. Thus, even for a large value of Vofc B , say 
0.70, the center of pressure moves only about 0.07 per- 
cent of the chord toward the origin. 
Again, if the angle of attack is increased to 4°, 

c c _ 1+0.54975 n 
c t 1+0.53656 11 

which shows that c e is nearer the nose of the profile 
than c { . The center of pressure in this case, with 
Vo/c a =0.70, moves about 0.13 percent of the chord 
toward the nose. 

This numerical example indicates that, although the 
effect of compressibility on the lift and the moment of 
a thin airfoil may be large, namely 



£ «-^=0.25517, Mc ~ M< =0.25847 



Mi 



for 



and 



0=0°, -=0.70 

Co 



U-L 



t 



\f Tut 

0.26938, ■ , r =0.26291 



Mt 



Ct 



l+|(M+2Vsin 2 j8) 



(32) 



Again, as in the case of elliptic cylinders, a neutral 
value for 0 is obtained when the centers of pressure in 
the compressible and the incompressible fluids coincide. 
It may be shown by means of numerical examples that, 
when the angle of attack is less than the neutral value 
of j8, the center of pressure c e in the compressible fluid 
moves from its position c t in the incompressible fluid 
toward the origin and that, when the angle of attack 
is greater than the neutral value of /S, the center of 
pressure c c moves toward the nose of the profile. 



for 



0=4°, ^=0.70 
Co 



the effect on the center of pressure may be considered 
negligible. 

As an example of a thicker airfoil, let e=0.10. The 
thickness of JoukowsH airfoils is proportional to e and 
a value of e of 0.10 gives a maximum thickness of about 
0.13 times the chord, a value rarely exceeded in practice. 
For this case, equation (32) becomes: 
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„ 1 +5(5.43168-4.34449 cos 2)3) 

c c_. & 

C< 1 +^(2.55446-1.44293 cos 2/8) 

Putting c c /c<=l yields 

cos 2)3=0.99161 

or a neutral angle 

/3=3°42' 

If the angle of attack is taken to be 0° and ^/co=0.70, 
then 

C 1.26636 - pfl 
c, 1.27232 -°- 99531 

or the center of pressure is displaced toward the origin 
a distance equal to 0.12 percent of the chord. 

If the angle of attack is increased to 6° with 
%/co=0.70, then 

Cc^ 1.28962 
Ci 1. 28005 : 



4.00748 



or the center of pressure is displaced toward the nose 
of the airfoil a distance equal to 0.19 percent of the 
chord. 

CONCLUDING REMARKS 

It has been shown that, for a thin airfoil and for 
small angles of attack, the moment is given by equation 
(31): 

M c =M«(l+|) 

This result is, however, limited to small values of (t. 
On the other hand, Glauert (reference 3) has shown 
that the lift on an airfoil in a compressible fluid is 
given by 

=£<(i+f+ • • •) 

The validity of this formula depends on the assumption 
that the velocity at the surface of the airfoil does 
not differ appreciably from the undisturbed velocity Vo, 
although the velocity in the undisturbed stream may 
be large. This assumption means that the airfoil must 
be thin and the angle of attack small. The similarity 
of the two preceding formulas, insofar as the first power 
of ii is concerned, strongly suggests that the expression 
for the moment on a thin airfoil at small angles of attack 
and for large stream velocities is given by: 



(33) 



This result is indeed implied in Glauert's work, where 
it is stated that the lift distribution along the chord 
remains unaltered but that the strength of each ele- 
mentary vortex is increased by the factor (1— ja)~*. 

In view of this discussion, it would appear that a 
more accurate expression for the moment M e may be 



obtained by introducing the factor (1— /t) -1 into equa- 
tion (20). Thus, equation (20) is replaced by 



\M t + poVoATm cos (0+ff) 



Vl-M 
+f[~§AA>F t m cos 03+O-) 



-2WR.P.^S^ 5 } 



(34) 



This formula differs from equation (20) in that it is 
valid for large stream velocities and differs from equa- 
tion (33) in that it estimates the effect of large dis- 
turbances to the main flow Vq because of the presence 
of an obstacle. 

As an example, consider the case of symmetrical 
Joukowski airfoils with sharp trailing edges. Equa- 
tion (26) is then replaced by 



Mc VT^( 1+ i+«(i+e) r, 

Mf 3 6(1+6) (l+e)» fT -1) . 

+!L~2 l + 6 (l +e ) + l + e(l + 6) (7+J) J (35 > 

In an analogous manner, the formula for the lift becomes 
L c =-JzLJj\+§(-l+M+N sin*/3)] (36) 

As a numerical example, consider an airfoil for which 

6= 0. 10. For this case 

M=1.11153 
#=2.88585 

Ar 



=i-=g(1.11153+2.8S5S5 sin'p) 



and 



(l+e) 8 (J+J) =4.80317-4.66366 cos 2j3 

2 Ar 

Table V gives the values of - p- and of (l+e)*(J+J) 

for several angles of attack. 

TABLE V 



p 


s Ar 


(1+«)«U+J) 


(deg.) 


»r f 




0 


1.11163 


0.13S81 


2 


1.11605 


.16089 


4 


1. 12667 


. 18488 


8 


1.14306 


.2IU1 


S 


1. 1674S 


.320)8 


10 


1.13S56 


.42078 



Table VI lists the values of the ratios M a /M, and 
L e [L t as given by equations (35) and (36), respectively, 
for various values of t>o/c 0 . The last column gives the 
values of MJMi (or LJLi) calculated from equation 
(33). 
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M. 

Mi 


L. 

T, 


1 




























Maw \ 

(deg.J 


0 


2 


i 


6 


8 


10 


0 


2 


4 


6 


8 


10 




0 

.40 
.50 
.80 
.70 
.7S 
.SO 
.85 


i.ooono 

1. 09871 
1.10732 
1.26362 
1.43005 
L 54898 
1.71317 
L 98812 


1.00000 
L 09964 
1.1HS85 
1.27200 
1.43383 
1.55347 
1.71S82 
I.9853S 


1.00000 
1.10240 
1. 17342 
1.27913 
1.44469 
L 66693 
L 73571 
1.98710 


1.00000 
1.10700 
1. 18102 
1.2909S 
1.46276 
1.58933 
1.76379 
2. 02821 


LOOOOO 
L11340 
1. 19168 
L 30749 
L48794 
L 62053 
L80293 
2.07353 


1.00000 
L12165 
1.20518 
L 32857 
1.82008 
L 66038 
1.85291 
2.13779 


1.00000 
1. 10084 
1. 17083 
1.27609 
1.43884 
L 86082 
L72616 
L 97480 


1.00000 
1. 10114 
1.17184 
1.276S9 
1.43975 
1.560*1 
L 72803 
L97722 


i.onooo 

L 10206 

l.iTzse 

I.27S25 
1 44336 
L 86829 
1.73364 
LSS443 


1.00000 
1. 10359 
LI753S 
1.28219 
1.44936 
1.57272 
L 74287 
1.99643 


1.00000 
1.10672 
1. 17890 
1.28767 
1.48772 
1.58303 
1.75696 
2.01314 


1.00000 
1.10843 
1.18339 
1.29467 
1.46840 
1.89682 
1.77266 
2.03448 


1.00000 
1.09110 
1. 16478 
1.26000 
1.40028 
L811S9 
1.66667 
I.S9S32 



Figures 5 and 6 show graphs based on table VI with 
MJMt and LJL ( , respectively, as ordinates and the 
Mach number Vo/ca as the abscissa. The dashed curves 
represent the Glauert approximation 

M e== L. = 1 

Mi r { VT^Ti 

It is seen from an examination of the table and the 
curves that, below the neutral angle of attack (in this 
case 3°42'), the Glauert approximation is better for the 
moment than for the lift but that, above this angle, the 
approximation is better for the lift than for the mo- 

2£ 
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Mi 
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'l.6 
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Fiouwt 6.— The variation of the ratio MJMi with the Maoh number itia. 
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Fioubi 6.— The variation of the ratio LJLt with the Mach number tt/et. 

ment. In any case, it appears that, at least for thin 
airfoils and for small angles of attack, the Glauert ap- 
proximation for both the moment and the lift is suffi- 
ciently good. It then follows that the effect of com- 
pressibility on the center of pressure is negligible. 
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APPENDIX A 



M= 



ds 



EFFECT OF COMPRESSIBILITY ON THE COUPLE ABOUT 
AN ELLIPTIC CYLINDER 

The moment on a body due to the fluid motion is, 
according to equation (1), given by 

f=- j[m]p 

where the positive direction of the normal vector n is 
taken from the body into the fluid. The components 
of the vectors r and n are, respectively, (x, y) and (cos nr, 
cos ny). 
Therefore 

[r n]=a; cos ny—y cos nx 
From figure 1, it is seen that 

dx=ds cos ny and dy=—ds cos nx 

Hence 

M=-^>p(x dx+y dy) = -\^p dr 2 

Now, the equations of transformation from Cartesian 
to elliptic coordinates are: 

x=c cosh £ cos rj, y=c sinh £ sin 17 

or 

s—c cosh (£+117) 

where £ takes on all values from zero to infinity and n 
takes on all values from 0 to 2t. Then £= constant 
and 77= constant represent confocal elhpses and hyper- 
bolas, respectively, the distance between the foci being 
2c. 

For any given ellipse £=a, say, 

x=c cosh a cos 17, y=c sinh a sin 17 

or 

r s =c 2 (cosh 2 a cos 2 q+sinh 2 a sin* ij) 

Then 

dr*=—c? sin 2rj dij 

Therefore 



1 f 2r 
M= — ^c 2 1 f sm2ij dr] 



Now, according to equation (7), the pressure p in a 
compressible fluid is given by 

p=constant ipot),, 2 — ^pof t 2 +^^-+ ... 
so that 

M c ==^poC 2 I v c 2 sin 2r t dri+yxp<fi i \ v, 2 sin 2^ 

The last two integrals in this expression are easily 
obtained since they involve a knowledge only of the 
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velocity » t in the incompressible fluid. Thus, suppose 
the elliptic cylinder to be in a flow of velocity t\> inclined 
at an angle 0 to the negative direction of the real axis 
and the circulation to be taken as zero. Then, if r„ is 
the radius of the circle into which the ellipse f=a is 
mapped by the Joukowski transformation, 

and the complex potential of the flow past the circle is 
it follows that the complex velocity is given by 



1 dw_ 
Vo~5z~ 



e^- 



n e vt 



z"e 



1- 



Atf* 



On the surface of the cylinder, z'=r^*. Therefore, 

Va dT . i 



4r 0 ; 



and 



where 



v?^? l-cos2(g+/S) 
Vo 2 "l-2«7 3 cos26i+(r 4 



"2r 0 



Now, from the Joukowski transformation, 



z'= 



where, the positive sign of the radical has been chosen 
in order that the regions at infinity of the 2 and s' 
planes shall coincide. 

Then if (j-, y) is replaced by the elliptic coordinates 
(£, rj), i. e., g=c cosh (£+117), it follows that: 



or, on the surface of the cylinder, 



Therefore 
and 



a=e~« and 6=y 



v? l-cos2(77+/3) y«[l-cos2(i;+)3)] 

1— 2e -2a cos 217+6"*" cosh 2a— cos 217 
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With this expression for VtjvJ, it is very easy to calcu- 
late the last two integrals in the expression for M e . 



Thus 



1 f ,r 

g/ijoocM t?, s sin 2,]dv 



Now 



and 



= 8^ W °Jo cosh2a-cos 2> 2 ^ 



X 



sin 27ii} drj 



X 



o cosh 2a— cos 2i?~ 



o cosh 2a— cos 2ij 



Therefore 



Similarly, 



1 fir i 

-gjupoC* I v? sin 2ij<f^=^jua-ptf»QV sin 2/3 



and 



X 



sin 2w?; drt 



o (cosh 2a— cos 2ij) s 



X 



cos 2nn d n 2ve _ itta y 2a ( th 2a+n) 
o (cosh 2 a— cos 2rj) 2 



Hence, 

■y*aflW e -'-_e» a 8111 2 " 
The sum of the two integrals becomes simply 

1 , .cosh 2a— cos 2/3 . na 
8111 2 ^ 

It is much more difficult to calculate the first iutegral 
in the expression for M e , for it is necessary in this case 
to know the velocity v c at the surface of the elliptic 
cylinder in a compressible fluid. For this purpose, it 
is convenient to make use of equation (13) of reference 
2. This equation, when applied to the case of an 
elliptic profile, becomes 

Ap jtfcos 2/3—0* . , . 

+S (2*1+1) sin [(2ji+1),,+ /S] fx ! »+WX 
«-o Jo 

-S(2n+1) sin [(2n+l)ij— /3] Px^afa+arfX 

n-0 JO 
» pi 

+4 sin /32> cos 2nv X^a^dX 

n-l Jo 

—sin 2/32ff s " -2 cos [(2n+l)v+0\ 

N-l 

+sin 2/3]Ty- J cos [(2n— /3] 



+2 sin 0 sin 20J2o 2n -* sin 2nr, 



where X=r 0 /r, c=c/2r a , and the a^'s are the coefficients 
of the cosine terms in the Fourier development of vf/itf 
in the region of flow. For the elliptic cylinder 

r t ii _ l-2X 2 cos 2(fl+/3)+X* 
l-2o 2 X I cos2e+(7*X* 



so that 



=o«o+S fe» cos 2ne+6i„ sin 2nd) 



_ n l + \ i Q.-2a t cos 2/3) 



o 0 =2 



1-<^X* 



n _„ g*(l+X 4 )-(l + q*X 4 ) cos 2)3 „ 
a*- 2, ^(l-^x*) ^ 

62,= 2 ^(ffX) 2a sin 2/3 

Now, it is recalled that Avfvo refers to the plane of the 
eircle into which the elliptic profile is mapped by the 
Joukowski transformation. Therefore 

or 

. „ sin 2t] r„ . . . „. , Av~[ 

i? Sm 2 ^ l-2a* C os2 v +a{ 2 Sm 

But 

— a**-* sin 2n v 

1— 2<r cos 2tj+tr* j=i ' 

Therefore, the first integral in the expression for M„ 
becomes 

iAftVj^ sin V sin 2^{l +f [ C0S 2 J~ 

+2 i-^cos2g ^ (27t _ r xs »-i ais _ 2£fx 

1 — ff* n=l Jo 

-i}(2n+ 1)0*-» J V^^]) 

The last two integrals in this expression may be written 
as 

2Snr»-« r^'-'idin-i-a^dX 
n-l Jo 

f V^-K^-j+a^dX 

n-l JO 

=22^" -6 fVa+^-a+^X 4 ) cos 2/S]X 4fl - J dX i 

n-l Jo 

_-fy»-6 r V(l + X«)-(l + a A X«)cos2g 

n-I Jo 1 — o^X 4 

_ 9 r Vq + X^-a + ^X*) cos 2/3 , 
~ 2 Jo ^(l-a'X 4 )* dX 

-X V( ' +x Vitff C03 ^ +^- 



-(l+ff^X^-W 



_ f V(l+X 4 )-(l+(7 4 X 4 ) cos2/3 ^ !! _ 1 
"Jo a*(I-<r*X*} tf X 



— a 2 cos 2/3 



1— 2«7 a cos 2/3+ff 4 , 1+a* 
2? log I=^ 
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where 



lfW<?f 0 *tf sin 2^=f p^V sin 2)3{l + |[ cos 2 jjT* 
-g*cos 2)3 



, - 1— a 3 cos 2/3 1— g*cos 2)3 



1-<T* 

1 -2(^003 2)3+ g* . 

2<7 9 



Filially 

. 1-2<t* cos 2jS+o-* 



pf cos 2)3— J 2 , n l— g* cos 2)3 1 — (r 2 cos 2)3 



2 ff 6 



1 +0* 1 cosh 2a— cos 2ft 



THE CALCULATION OF Ao,/po FOR AN ELLIPTIC 
CYLINDER AT 

Equation (13) of reference 2 gives an expression for 
Avi/va independent of the shape of the obstacle. In 
order to evaluate A»j/» 0 , it is necessary to know the 
Fourier development of VijvJ. For an elliptic cylinder, 
with the circulation chosen so that the rear stagnation 
point occurs at the end of the major axis, it may be 
easily shown by means of the Joukowski transformation 
and the complex velocity potential of the flow about a 
circular cylinder that: 

v t * (1+2X cos d+\*)[l— 2X cos (g+2p)+X*I 
n? l-2o 3 X ! cos 26+<r*\* 

1 °° 

== s ee o+S («» cos sin 7i0) 

where 

X=r 0 /r, <r=c/2rQ, )3 is the angle of attack, and r, 6 
are the polar coordinates of a point in the region of 
flow. 

The Fourier series for Cj 1 /^ 2 is obtained by making use 
of the expansion 



1 



^=r=W<[i+2g(^rcos2^] 



l-2o 2 X s cos20+ff' 

and it may be shown without difficulty that 



_ 2(1 + X 2 ) a 4X* cos 2)3 
a "~ l-o*X* 1-^X 1 



l + ^X 1 



^(l-^X 8 ) 

_ 4X(1 + X') sin*/3 , 
Osn+i— 1 _ g a x a (<rX) 2 " 

6ii»=2X J (o-X) s *- 2 sin 2/3 

, _ 2XQ + X*) sin 2)3 
RVX 1 



cos 2/3 



(<A) S 



A^*=f Po«o s c* sin 2/3 
Replacing e~ Sa by a 2 : 

+f (l-2o* cos 20+^ log J,)] 

This expression agrees with equation (21) when the 

sinh 2a j circulation is taken to be zero. 

APPENDIX B 

Equation (13) of reference 2 may then be written as 

\ » 0 /j_ T „_i jo 

-a cos fx^^+^x 

B-U JO 

-cos 0Sf2n+l) fV^-^+arfX 

n=0 Jo 

+2 cos f'x^+adX-sin ,8 PxctotfX 

n-o Jo Jo 

-sinj9S(2ti+l) fV^+WX 

n-l Jo 

+2 sin f 1 X 2 <»+i>a 2 „ +1 <fX 
n-o Jo 

-sin )8SC2»+1) f'x^-W+s^X 

n-0 JO 

+2 sin fx^+^X 

n-0 Jo 

+4-sin |8ZJ(»+1) f'x^B+srfX 

n=0 Jo 

-2 sin 0S(2n+l) f "x^+idxl 

n-0 Jo J 

If the expressions for a n and are inserted and the 
integrations performed, it follows after considerable 
but straightforward labor that 

(~) s _=^sm f}-^ sin 3/?+^ [(l+«r)»Iog (l+«r) 

— (1— <^) a log (l-«r)] (sin /3+sin 3)3) 

-74 I(1+O j log (l+^-d-ff 3 )* log (l-a^sin jB 

+ 8l °gg- qA) sin')3-^sin)3log (l+o 2 ) 



](<rX)*» 



A THE ORETICAL STUDY OF THE MOMENT ON A BODY IN A COMPRESSIBLE FLUID 

The additional circulation AK is given by 
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(See section on the elliptic profile.) 
Therefore, with K f =2 sin P, it follows that 



where 



M=^+^ log (l + O 



+ —j? — lo s r=? — ?~ log w 



l + <r 



^=i±^logi±f-^[l+log (l+O] 



It is interesting to note that, for the limiting case of 
a circular cylinder <r=0, the foregoing equation yields 

and compares with Poggi's result (reference 5) 
/ AT\ = 11 

For the straight-line profile <7=1, it is seen that M=l 
and 2V=- 00 . For an infinitely thin profile and a vanish- 
ingly small angle of attack, however, 

and compares with Glauert's formula (reference 3): 



AT 

r, 



= — 1 + 



EVALUATION OF / 



APPENDIX C 

Therefore 



2n +l A„A a +i _ -y^ A n A a +i 



ta-l 



Since 

=log (!-#), ^!=2ir 0 sin ft At^aW-rtfer* 

n-l W 

and A t =a i ra( L ke^—e~^) l it follows that 

= ~$ HST^ 1-2 * cos 2/? + P ) lo = 
-^(l-2i cos 2/S+* s )(l +§#+§£*) 

-3i sin 



1— & , a 1+A 

e== r+& and r 0 = ^- 



Therefore 



7=R. P. 



-e^ 2fl+l ^l a+1 



2 sin 2/3^in(w+l) r 0 2n+1 
^p(l-2* cos 2/3+F)(l+|p+^) 



2n 



+ 
Now 



^(9+5A)(l-cos 2^) +1 4 



209142—40 



a 1 , , 1— 6 

r n 1+e 1 + e 

-so 



Hence 

(l+e)*/=£^?(l-2* cos 2/S+**)[^log ^ 
-(l+|t«+|i*)]+J(9+tt)(l-008 20) 

+^(3+£)(l-A) s 
For small values of or thick airfoils, 

(l+e)*/=(l+*) 2 (l-2A cos 2ft+t*)fyi+£tf* 

+ ik kl+ • ' •)+£p+5*)Cl-eos20 

+^(3+A)(l-^) ! 

EVALUATION OF J 

A n .i4 a+ 2 AiAj AgAj ■^- y A H -iA %+ i 

^i(n+l)r^ 2r 0 * + 3r 0 ° 7ir 0 2 » 

a«(l-2fc cos ^fc 2 " .a 2 . _„ „ 

= "Sir- v 3111 m 6 } 

- gt(1 - 2 ^ 2W) ( 1+ f4>v sin ^"-^ 
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Therefore 



J=-R. P. 



A. n A.- 



sin2^(ji+l)r 0 Sa+2 

g 4 (l-24 cob 20+4*) /, 4 s 4«\ 
+^- 2 (3+24)(l-24 cos 2/3+4*) 

+ g* (1 -^( r f*) 



and 

(l+e)V=£^(l-2* cos 2/S+4»)[| 5 log j^p 

-(l+|-+^)]+i(3+24)(l-24 cos 2/S+4 3 ) 

+A4(l-4)+^4 i (l-4)* 
For small values of 4, or thick airfoils, 
(l+e) 2 J=(l+4)*(l-24 cos 2/S+fc»)Qp 

.)+J(3+24)(l-24 cos 2/3+4*) 



+ 20 i4+ 24 i8+ 
+ i*(l-fc)+l*»(l. 



APPENDIX D 



EVALUATION OF (ABj/DoJa-r FOR THE CASE OF A 
JOUKOWSKI PROFILE 

According to equation (22) of reference 2, (A» 4 /ob)j_ 
takes the following form: 

+^^*(-3-34A+|4-A 2 4 2 +|4A 2 ) 
+ ?A^|- 1 _l iI+ I^ log(1 _ H) ] 

+ (1-A)*+2A(1-A) log (1+4)] 
+ 4A»(l-A)* pog Cl _ A)+4A »i og 



4A : 
4 

8(1 



1+7 



2(1-^) 

" 4* 



log (1-4 2 )] 



F ^- 3 [M 2 - 1 -2A f( i_ A ) 2+ 2A(i- A ) log (1+4)] 



-^p^-Vg (l-4)+4A 2 log (1+4) -4 s ] 



where 



A=r4— and 4=? * 



"l+e" 1 + e 

When this expression for (A^/co)*-* was obtained, a 
slight error was found in equations (19) of reference 2. 
The expressions for a 1 —a l and h—Hi should be as 
follows: 



-_ 2\ cos 2/S| 

d\ — £li £2" 



! [a(A+24)+2A»4X 2 ] 



and 



2A*\ I + X*(l-24cos2/3) 



61-61 2X ^ 2/? [a(A+24) +24A s X»] 



It is to be noted in tbe expression for (At>j/t%)i,., that 
most of the terms contain powers of 4 in the denomi- 
nator. It appears at first, then, that the coefficients 
of sin 18 and sin 8 /3 may become infinite for 4=0. This 
apparent difficulty disappears, however, when (Ai>j/t*))i-* 
is expressed as a power series in 4. It is then found 
that the terms involving reciprocal powers of 4 cancel 
and the following expression results: 



(-) 



=% sin 



309, 289,, 7 j, 
"160 672 + 32^" 



191 
'640' 



0 

2 8111 \ 24 240* 24^ 1680* 560* " " J 
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